Ancestral aneuploidy and stable chromosomal duplication resulting in differential genome structure and gene expression control in trypanosomatid parasites.

Aneuploidy is widely observed in both unicellular and multicellular eukaryotes, usually associated with adaptation to stress conditions. Chromosomal duplication stability is a tradeoff between the fitness cost of having unbalanced gene copies and the potential fitness gained from increased dosage of specific advantageous genes. Trypanosomatids, a family of protozoans that include species that cause neglected tropical diseases, are a relevant group to study aneuploidies. Their life cycle has several stressors that could select for different patterns of chromosomal duplications and/or losses, and their nearly universal use of polycistronic transcription increases their reliance on gene expansion/contraction, as well as post-transcriptional control as mechanisms for gene expression regulation. By evaluating the data from 866 isolates covering seven trypanosomatid genera, we have revealed that aneuploidy tolerance is an ancestral characteristic of trypanosomatids but has a reduced occurrence in a specific monophyletic clade that has undergone large genomic reorganization and chromosomal fusions. We have also identified an ancient chromosomal duplication that was maintained across these parasite's speciation, named collectively as the trypanosomatid ancestral supernumerary chromosome (TASC). TASC has most genes in the same coding strand, is expressed as a disomic chromosome (even having four copies), and has increased potential for functional variation, but it purges highly deleterious mutations more efficiently than other chromosomes. The evidence of stringent control over gene expression in this chromosome suggests that these parasites have adapted to mitigate the fitness cost associated with this ancient chromosomal duplication.

RAD51-mediated R-loop formation acts to repair transcription-associated DNA breaks driving antigenic variation in Trypanosoma brucei.

RNA-DNA hybrids are epigenetic features of all genomes that intersect with many processes, including transcription, telomere homeostasis, and centromere function. Increasing evidence suggests that RNA-DNA hybrids can provide two conflicting roles in the maintenance and transmission of genomes: They can be the triggers of DNA damage, leading to genome change, or can aid the DNA repair processes needed to respond to DNA lesions. Evasion of host immunity by African trypanosomes, such as Trypanosoma brucei, relies on targeted recombination of silent Variant Surface Glycoprotein (VSG) genes into a specialized telomeric locus that directs transcription of just one VSG from thousands. How such VSG recombination is targeted and initiated is unclear. Here, we show that a key enzyme of T. brucei homologous recombination, RAD51, interacts with RNA-DNA hybrids. In addition, we show that RNA-DNA hybrids display a genome-wide colocalization with DNA breaks and that this relationship is impaired by mutation of RAD51. Finally, we show that RAD51 acts to repair highly abundant, localised DNA breaks at the single transcribed VSG and that mutation of RAD51 alters RNA-DNA hybrid abundance at 70 bp repeats both around the transcribed VSG and across the silent VSG archive. This work reveals a widespread, generalised role for RNA-DNA hybrids in directing RAD51 activity during recombination and uncovers a specialised application of this interplay during targeted DNA break repair needed for the critical T. brucei immune evasion reaction of antigenic variation.

The dynamic subcellular localisation of Rad1 is cell cycle dependent in Leishmania major.

The subcellular localisation of Rad1, a subunit of the Leishmania major 9-1-1 complex, remains unexplored. Herein, we reveal that Rad1 localises predominantly to the nucleus. Upon hydroxyurea treatment, the diffuse nuclear localisation of Rad1 becomes more punctate, suggesting that Rad1 is responsive to replication stress. Moreover, Rad1 localisation correlates with cell cycle progression. In the majority of G1 to early S-phase cells, Rad1 localises predominantly to the nucleus. As cells progress from late-S phase to mitosis, Rad1 relocalizes to both the nucleus and the cytoplasm in ∼90 % of cells. This pattern of distribution is different from Rad9 and Hus1, which remain nuclear throughout the cell cycle, suggesting Leishmania Rad1 may regulate 9-1-1 activities and/or perform relevant functions outside the 9-1-1 complex.

Immunoprecipitation of RNA-DNA hybrid interacting proteins in Trypanosoma brucei reveals conserved and novel activities, including in the control of surface antigen expression needed for immune evasion by antigenic variation.

RNA-DNA hybrids are epigenetic features of genomes that provide a diverse and growing range of activities. Understanding of these functions has been informed by characterising the proteins that interact with the hybrids, but all such analyses have so far focused on mammals, meaning it is unclear if a similar spectrum of RNA-DNA hybrid interactors is found in other eukaryotes. The African trypanosome is a single-cell eukaryotic parasite of the Discoba grouping and displays substantial divergence in several aspects of core biology from its mammalian host. Here, we show that DNA-RNA hybrid immunoprecipitation coupled with mass spectrometry recovers 602 putative interactors in T. brucei mammal- and insect-infective cells, some providing activities also found in mammals and some lineage-specific. We demonstrate that loss of three factors, two putative helicases and a RAD51 paralogue, alters T. brucei nuclear RNA-DNA hybrid and DNA damage levels. Moreover, loss of each factor affects the operation of the parasite immune survival mechanism of antigenic variation. Thus, our work reveals the broad range of activities contributed by RNA-DNA hybrids to T. brucei biology, including new functions in host immune evasion as well as activities likely fundamental to eukaryotic genome function.

AAK1-like: A putative pseudokinase with potential roles in cargo uptake in bloodstream form Trypanosoma brucei parasites.

Selection and internalization of cargo via clathrin-mediated endocytosis requires adaptor protein complexes. One complex, AP-2, acts during cargo selection at the plasma membrane. African trypanosomes lack all components of the AP-2 complex, except for a recently identified orthologue of the AP-2-associated protein kinase 1, AAK1. In characterized eukaryotes, AAK1 phosphorylates the µ2 subunit of the AP-2 complex to enhance cargo recognition and uptake into clathrin-coated vesicles. Here, we show that kinetoplastids encode not one, but two AAK1 orthologues: one (AAK1L2) is absent from salivarian trypanosomes, while the other (AAK1L1) lacks important kinase-specific residues in a range of trypanosomes. These AAK1L1 and AAK1L2 novelties reinforce suggestions of functional divergence in endocytic uptake within salivarian trypanosomes. Despite this, we show that AAK1L1 null mutant Trypanosoma brucei, while viable, display slowed proliferation, morphological abnormalities including swelling of the flagellar pocket, and altered cargo uptake. In summary, our data suggest an unconventional role for a putative pseudokinase during endocytosis and/or vesicular trafficking in T. brucei, independent of AP-2.

The dynamic subcellular localisation of Rad1 is cell cycle dependent in Leishmania major

The subcellular localisation of Rad1, a subunit of the Leishmania major 9-1-1 complex, remains unexplored. Herein, we reveal that Rad1 localises predominantly to the nucleus. Upon hydroxyurea treatment, the diffuse nuclear localisation of Rad1 becomes more punctate, suggesting that Rad1 is responsive to replication stress. Moreover, Rad1 localisation correlates with cell cycle progression. In the majority of G1 to early S-phase cells, Rad1 localises predominantly to the nucleus. As cells progress from late-S phase to mitosis, Rad1 relocalizes to both the nucleus and the cytoplasm in ∼90 % of cells. This pattern of distribution is different from Rad9 and Hus1, which remain nuclear throughout the cell cycle, suggesting Leishmania Rad1 may regulate 9-1-1 activities and/or perform relevant functions outside the 9-1-1 complex.

Emergence and adaptation of the cellular machinery directing antigenic variation in the African trypanosome.

Survival of the African trypanosome within its mammalian hosts, and hence transmission between hosts, relies upon antigenic variation, where stochastic changes in the composition of their protective variant-surface glycoprotein (VSG) coat thwart effective removal of the pathogen by adaptive immunity. Antigenic variation has evolved remarkable mechanistic complexity in Trypanosoma brucei, with switching of the VSG coat executed by either transcriptional or recombination reactions. In the former, a single T. brucei cell selectively transcribes one telomeric VSG transcription site, termed the expression site (ES), from a pool of around 15. Silencing of the active ES and activation of one previously silent ES can lead to a co-ordinated VSG coat switch. Outside the ESs, the T. brucei genome contains an enormous archive of silent VSG genes and pseudogenes, which can be recombined into the ES to execute a coat switch. Most such recombination involves gene conversion, including copying of a complete VSG and more complex reactions where novel 'mosaic' VSGs are formed as patchworks of sequences from several silent (pseudo)genes. Understanding of the cellular machinery that directs transcriptional and recombination VSG switching is growing rapidly and the emerging picture is of the use of proteins, complexes and pathways that are not limited to trypanosomes, but are shared across the wider grouping of kinetoplastids and beyond, suggesting co-option of widely used, core cellular reactions. We will review what is known about the machinery of antigenic variation and discuss if there remains the possibility of trypanosome adaptations, or even trypanosome-specific machineries, that might offer opportunities to impair this crucial parasite-survival process.

Single-cell transcriptomic analysis of bloodstream Trypanosoma brucei reconstructs cell cycle progression and developmental quorum sensing.

Developmental steps in the trypanosome life-cycle involve transition between replicative and non-replicative forms specialised for survival in, and transmission between, mammalian and tsetse fly hosts. Here, using oligopeptide-induced differentiation in vitro, we model the progressive development of replicative 'slender' to transmissible 'stumpy' bloodstream form Trypanosoma brucei and capture the transcriptomes of 8,599 parasites using single cell transcriptomics (scRNA-seq). Using this framework, we detail the relative order of biological events during asynchronous development, profile dynamic gene expression patterns and identify putative regulators. We additionally map the cell cycle of proliferating parasites and position stumpy cell-cycle exit at early G1 before progression to a distinct G0 state. A null mutant for one transiently elevated developmental regulator, ZC3H20 is further analysed by scRNA-seq, identifying its point of failure in the developmental atlas. This approach provides a paradigm for the dissection of differentiation events in parasites, relevant to diverse transitions in pathogen biology.

Unpicking the Roles of DNA Damage Protein Kinases in Trypanosomatids.

To preserve genome integrity when faced with DNA lesions, cells activate and coordinate a multitude of DNA repair pathways to ensure timely error correction or tolerance, collectively called the DNA damage response (DDR). These interconnecting damage response pathways are molecular signal relays, with protein kinases (PKs) at the pinnacle. Focused efforts in model eukaryotes have revealed intricate aspects of DNA repair PK function, including how they direct DDR pathways and how repair reactions connect to wider cellular processes, including DNA replication and transcription. The Kinetoplastidae, including many parasites like Trypanosoma spp. and Leishmania spp. (causative agents of debilitating, neglected tropical infections), exhibit peculiarities in several core biological processes, including the predominance of multigenic transcription and the streamlining or repurposing of DNA repair pathways, such as the loss of non-homologous end joining and novel operation of nucleotide excision repair (NER). Very recent studies have implicated ATR and ATM kinases in the DDR of kinetoplastid parasites, whereas DNA-dependent protein kinase (DNA-PKcs) displays uncertain conservation, questioning what functions it fulfills. The wide range of genetic manipulation approaches in these organisms presents an opportunity to investigate DNA repair kinase roles in kinetoplastids and to ask if further kinases are involved. Furthermore, the availability of kinase inhibitory compounds, targeting numerous eukaryotic PKs, could allow us to test the suitability of DNA repair PKs as novel chemotherapeutic targets. Here, we will review recent advances in the study of trypanosomatid DNA repair kinases.

Application of single-cell transcriptomics to kinetoplastid research.

Kinetoplastid parasites are responsible for both human and animal diseases across the globe where they have a great impact on health and economic well-being. Many species and life cycle stages are difficult to study due to limitations in isolation and culture, as well as to their existence as heterogeneous populations in hosts and vectors. Single-cell transcriptomics (scRNA-seq) has the capacity to overcome many of these difficulties, and can be leveraged to disentangle heterogeneous populations, highlight genes crucial for propagation through the life cycle, and enable detailed analysis of host­parasite interactions. Here, we provide a review of studies that have applied scRNA-seq to protozoan parasites so far. In addition, we provide an overview of sample preparation and technology choice considerations when planning scRNA-seq experiments, as well as challenges faced when analysing the large amounts of data generated. Finally, we highlight areas of kinetoplastid research that could benefit from scRNA-seq technologies.

The MRN complex promotes DNA repair by homologous recombination and restrains antigenic variation in African trypanosomes.

Homologous recombination dominates as the major form of DNA repair in Trypanosoma brucei, and is especially important for recombination of the subtelomeric variant surface glycoprotein during antigenic variation. RAD50, a component of the MRN complex (MRE11, RAD50, NBS1), is central to homologous recombination through facilitating resection and governing the DNA damage response. The function of RAD50 in trypanosomes is untested. Here we report that RAD50 and MRE11 are required for RAD51-dependent homologous recombination and phosphorylation of histone H2A following a DNA double strand break (DSB), but neither MRE11 nor RAD50 substantially influence DSB resection at a chromosome-internal locus. In addition, we reveal intrinsic separation-of-function between T. brucei RAD50 and MRE11, with only RAD50 suppressing DSB repair using donors with short stretches of homology at a subtelomeric locus, and only MRE11 directing DSB resection at the same locus. Finally, we show that loss of either MRE11 or RAD50 causes a greater diversity of expressed VSG variants following DSB repair. We conclude that MRN promotes stringent homologous recombination at subtelomeric loci and restrains antigenic variation.

Read, Write, Adapt: Challenges and Opportunities during Kinetoplastid Genome Replication.

The genomes of all organisms are read throughout their growth and development, generating new copies during cell division and encoding the cellular activities dictated by the genome's content. However, genomes are not invariant information stores but are purposefully altered in minor and major ways, adapting cellular behaviour and driving evolution. Kinetoplastids are eukaryotic microbes that display a wide range of such read-write genome activities, in many cases affecting critical aspects of their biology, such as host adaptation. Here we discuss the range of read-write genome changes found in two well-studied kinetoplastid parasites, Trypanosoma brucei and Leishmania, focusing on recent work that suggests such adaptive genome variation is linked to novel strategies the parasites use to replicate their unconventional genomes.

Genome duplication in Leishmania major relies on persistent subtelomeric DNA replication.

DNA replication is needed to duplicate a cell's genome in S phase and segregate it during cell division. Previous work in Leishmania detected DNA replication initiation at just a single region in each chromosome, an organisation predicted to be insufficient for complete genome duplication within S phase. Here, we show that acetylated histone H3 (AcH3), base J and a kinetochore factor co-localise in each chromosome at only a single locus, which corresponds with previously mapped DNA replication initiation regions and is demarcated by localised G/T skew and G4 patterns. In addition, we describe previously undetected subtelomeric DNA replication in G2/M and G1-phase-enriched cells. Finally, we show that subtelomeric DNA replication, unlike chromosome-internal DNA replication, is sensitive to hydroxyurea and dependent on 9-1-1 activity. These findings indicate that Leishmania's genome duplication programme employs subtelomeric DNA replication initiation, possibly extending beyond S phase, to support predominantly chromosome-internal DNA replication initiation within S phase.

Genome maintenance functions of a putative Trypanosoma brucei translesion DNA polymerase include telomere association and a role in antigenic variation.

Maintenance of genome integrity is critical to guarantee transfer of an intact genome from parent to offspring during cell division. DNA polymerases (Pols) provide roles in both replication of the genome and the repair of a wide range of lesions. Amongst replicative DNA Pols, translesion DNA Pols play a particular role: replication to bypass DNA damage. All cells express a range of translesion Pols, but little work has examined their function in parasites, including whether the enzymes might contribute to host-parasite interactions. Here, we describe a dual function of one putative translesion Pol in African trypanosomes, which we now name TbPolIE. Previously, we demonstrated that TbPolIE is associated with telomeric sequences and here we show that RNAi-mediated depletion of TbPolIE transcripts results in slowed growth, altered DNA content, changes in cell morphology, and increased sensitivity to DNA damaging agents. We also show that TbPolIE displays pronounced localization at the nuclear periphery, and that its depletion leads to chromosome segregation defects and increased levels of endogenous DNA damage. Finally, we demonstrate that TbPolIE depletion leads to deregulation of telomeric variant surface glycoprotein genes, linking the function of this putative translesion DNA polymerase to host immune evasion by antigenic variation.

Conditional knockout of RAD51-related genes in Leishmania major reveals a critical role for homologous recombination during genome replication.

Homologous recombination (HR) has an intimate relationship with genome replication, both during repair of DNA lesions that might prevent DNA synthesis and in tackling stalls to the replication fork. Recent studies led us to ask if HR might have a more central role in replicating the genome of Leishmania, a eukaryotic parasite. Conflicting evidence has emerged regarding whether or not HR genes are essential, and genome-wide mapping has provided evidence for an unorthodox organisation of DNA replication initiation sites, termed origins. To answer this question, we have employed a combined CRISPR/Cas9 and DiCre approach to rapidly generate and assess the effect of conditional ablation of RAD51 and three RAD51-related proteins in Leishmania major. Using this approach, we demonstrate that loss of any of these HR factors is not immediately lethal but in each case growth slows with time and leads to DNA damage and accumulation of cells with aberrant DNA content. Despite these similarities, we show that only loss of RAD51 or RAD51-3 impairs DNA synthesis and causes elevated levels of genome-wide mutation. Furthermore, we show that these two HR factors act in distinct ways, since ablation of RAD51, but not RAD51-3, has a profound effect on DNA replication, causing loss of initiation at the major origins and increased DNA synthesis at subtelomeres. Our work clarifies questions regarding the importance of HR to survival of Leishmania and reveals an unanticipated, central role for RAD51 in the programme of genome replication in a microbial eukaryote.

Targeting the trypanosome kinetochore with CLK1 protein kinase inhibitors.

The kinetochore is a macromolecular structure that assembles on the centromeres of chromosomes and provides the major attachment point for spindle microtubules during mitosis. In Trypanosoma brucei, the proteins that make up the kinetochore are highly divergent; the inner kinetochore comprises at least 20 distinct and essential proteins (KKT1-20) that include four protein kinases-CLK1 (also known as KKT10), CLK2 (also known as KKT19), KKT2 and KKT3. Here, we report the identification and characterization of the amidobenzimidazoles (AB) protein kinase inhibitors that show nanomolar potency against T. brucei bloodstream forms, Leishmania and Trypanosoma cruzi. We performed target deconvolution analysis using a selection of 29 T. brucei mutants that overexpress known essential protein kinases, and identified CLK1 as a primary target. Biochemical studies and the co-crystal structure of CLK1 in complex with AB1 show that the irreversible competitive inhibition of CLK1 is dependent on a Michael acceptor forming an irreversible bond with Cys 215 in the ATP-binding pocket, a residue that is not present in human CLK1, thereby providing selectivity. Chemical inhibition of CLK1 impairs inner kinetochore recruitment and compromises cell-cycle progression, leading to cell death. This research highlights a unique drug target for trypanosomatid parasitic protozoa and a new chemical tool for investigating the function of their divergent kinetochores.

Replication origin location might contribute to genetic variability in Trypanosoma cruzi.

BACKGROUND: DNA replication in trypanosomatids operates in a uniquely challenging environment, since most of their genomes are constitutively transcribed. Trypanosoma cruzi, the etiological agent of Chagas disease, presents high variability in both chromosomes size and copy number among strains, though the underlying mechanisms are unknown. RESULTS: Here we have mapped sites of DNA replication initiation across the T. cruzi genome using Marker Frequency Analysis, which has previously only been deployed in two related trypanosomatids. The putative origins identified in T. cruzi show a notable enrichment of GC content, a preferential position at subtelomeric regions, coinciding with genes transcribed towards the telomeres, and a pronounced enrichment within coding DNA sequences, most notably in genes from the Dispersed Gene Family 1 (DGF-1). CONCLUSIONS: These findings suggest a scenario where collisions between DNA replication and transcription are frequent, leading to increased genetic variability, as seen by the increase SNP levels at chromosome subtelomeres and in DGF-1 genes containing putative origins.

Trypanosoma brucei and Trypanosoma cruzi DNA Mismatch Repair Proteins Act Differently in the Response to DNA Damage Caused by Oxidative Stress.

MSH2, associated with MSH3 or MSH6, is a central component of the eukaryotic DNA Mismatch Repair (MMR) pathway responsible for the recognition and correction of base mismatches that occur during DNA replication and recombination. Previous studies have shown that MSH2 plays an additional DNA repair role in response to oxidative damage in Trypanosoma cruzi and Trypanosoma brucei. By performing co-immunoprecipitation followed by mass spectrometry with parasites expressing tagged proteins, we confirmed that the parasites' MSH2 forms complexes with MSH3 and MSH6. To investigate the involvement of these two other MMR components in the oxidative stress response, we generated knockout mutants of MSH6 and MSH3 in T. brucei bloodstream forms and MSH6 mutants in T. cruzi epimastigotes. Differently from the phenotype observed with T. cruzi MSH2 knockout epimastigotes, loss of one or two alleles of T. cruzi msh6 resulted in increased susceptibility to H2O2 exposure, besides impaired MMR. In contrast, T. brucei msh6 or msh3 null mutants displayed increased tolerance to MNNG treatment, indicating that MMR is affected, but no difference in the response to H2O2 treatment when compared to wild type cells. Taken together, our results suggest that, while T. cruzi MSH6 and MSH2 are involved with the oxidative stress response in addition to their role as components of the MMR, the DNA repair pathway that deals with oxidative stress damage operates differently in T. brucei.

Next-Generation Analysis of Trypanosomatid Genome Stability and Instability.

Understanding the rate and patterns of genome variation is becoming ever more amenable to whole-genome analysis through advances in DNA sequencing, which may, at least in some circumstances, have supplanted more localized analyses by cellular and genetic approaches. Whole-genome analyses can utilize both short- and long-read sequence technologies. Here we describe how sequence generated by these approaches has been used in trypanosomatids to examine DNA replication dynamics, the accumulation of modified histone H2A due to genome damage, and evaluation of genome variation, focusing on ploidy change.